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Abstract 



We show the existence of bound two-electron states in an almost depleted 
two-dimensional island. These two-electron states are carried by special com- 
pact configurations of four single-electron levels. The existence of these states 
does not require phonon mediation, and is facilitated by the disorder-induced 
potential relief and by the electron-electron repulsion only. The density of 
two-electron states is estimated and their evolution with the magnetic field is 
discussed. 
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In a recent experimentcra the tunneling of electrons from a semimetallic electrode into 
the localized states (LS's) in a quantum well was studied. By using the electron beam 
lithography, a dot with a diameter as asmall as 1 /iin was formed. The small size of the 
dot enabled the authors to detect the individual tunneling acts. These acts manifested 
themselves as narrow peaks in the differential capacitance measured as a function of the 
bias applied across the structure. To introduce the LS's in the GaAs quantum well the 
neighboring AlGaAs region was doped with Si donors. 

The authors identified the origin of LS's by studying the evolution of the peak positions 
with magnetic field, B, perpendicular to the well. They assocaited the peaks with electronic 
states of essentially two types: (i) ground state in the cylindrically symmetric parabolic 
potential; the .£>- dependence of the corresponding peaks above certain B approached the 
one for the lowest Landau level. Such a parabolic confinement presumably results from the 
fluctuations in the concentration of donors in the barrier, (ii) bound state of an electron 
at a Si donor which could migrate into the well during the growth. For these states the 
increase of the energy with B was much slower than for the group (i). 

There is a puzzling feature in the data reported in0: few of the peaks observed were 
twice as high as isolated one-electron peaks. This suggests that two electrons tunnel into 
the well at the same voltage applied. By tracing the -B-dependencies of these peaks the 
authors have ruled out the possibility that they reflect accidental degeneracies in the energy 
positions of LS's in some distant minima. All double peaks retained their height within a 
certain range of B and then split into doublets at some critical value of the magnetic field. 
Such a behavior indicates that both LS's involved "feel" each other and, thus, are located 
close in space. On the other hand it is apparent that two close potential minima cannot 
accomodate two electrons at the same bias. Even if the energy levels are degenerate, tun- 
neling of one electron would elevate the level for the second electron, so that the subsequent 
tunneling will occur at the bias larger by the energy of the Coulomb repulsion. The authors 
mentioned that the physical mechanism which could resolve this paradox is the polaronic 
effect: this effect favors double occupation of LS's in glassy materials!. Their conjecture 
was further developed in@, where the two-electron state in a hybrid hydrogenic-parabolic 
potential in the presence of electron-phonon interaction was considered. The pair-binding 
condition used inu implicitly assumed that the two electrons share the same lattice deforma- 
tion, which leads to the enhancement of the polaronic effect. In fact, the used in@ polaronic 
shift per electron in the paired state is twice the shift for a single localized electron. Under 
this assumption bound two-electron states were found even at weak electron-phonon inter- 
action, provided that the distance between the hydrogen-like impurity and the center of the 
parabolic potential is larger than 8a , where ao is the radius of the hydrogenic state. For 
such distances the enhanced polaronic shift overweighs the Coulomb repulsion. However we 
find the underlying assumption hard to justify. Indeed, the spatial scale of the polaronic 
deformation coincides with the size of a single-electron state!; two distant electrons do not 
share the same deformation, and, therefore, the corresponding enhancement of the polaron 
shift is suppressed. 

In the present paper we demonstrate that the double peaks observed in0 can be naturally 
explained without invoking the electron-phonon interactions. Our explanation is based ex- 
clusively on electrostatics. We assume that the electrons are strongly localized and neglect 
the overlap of their wave functions and, correspondingly, the exchange interaction. On the 
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other hand the modification of the Coulomb interaction between localized electrons due to 
the presence of an electrode plays an important role in our picture. We show that, in con- 
trast toi, the objects responsible for two-electron tunneling are the compact groups of LS's 
and calculate the relative portion of the double peaks. 

First of all, let us establish the general criterion for two-electron tunneling. Consider a 
cluster of N LS's occupied by n electrons. The distribution of electrons over LS's corresponds 
to the minimal possible energy which we denote as E^. The position of the Fermi level in the 
electrode, Ep, at which an additional electron will enter the cluster is determined from the 
condition Ep + E 1 ^ = E^ +1 . If two electrons enter the cluster, the corresponding position 
of the Fermi level, Ep, satisfies the relation 2Ep + E^ = E^ +2 . Double peak occurs if 
Ep < Ep. This leads us to the following criterion: 

E n N +2 + E n N <2E n N +l . (1) 

Obviously condition (Jl|) cannot be satisfied if AT = 2. Formally, if we denote the energies 
of two LS's as Si and e 2 so that E\ < e 2 , then E\ = e\, E\ = e% + e 2 + V\ 2 , V\ 2 being the 
interaction energy of two electrons occupying the first and the second LS's. We see that 
E\ > 2i?2-the inequality opposite to (|l|). 

Now we will prove that for N = 3 the occurrence of a double peak is also forbidden 
by Eq. ([[]). By analogy to the consideration above, the case n = is obvious for any 
N since a single electron on a cluster will occupy the LS with the lowest energy level, so 
that the condition ([I]) is violated even without the Coulomb repulsion. Thus the only case 
to be considered is n = 1. Let us again order the energies of LS's: £\ < e 2 < e 3 . Then 
£"3 = Si, and El = e\ + e 2 + £3 + V\ 2 + V13 + V 23 . Now there are three candidates for E 3 in 
accordance with three variants of occupation of iV = 3 cluster by two electrons. Important 
is that E\ is the minimal of these three energies. This means that if the double peak is 
possible, the condition (|l|) should be met when we substitute for E\ each of these energies. 
Let us choose two of the three candidates for E\, namely: e\ + e 2 + V\ 2 and e% + £3 + V13, 
which correspond to the occupation of the first and second, and the first and the third LS's 
respectively. Substituting them into the right-hand side of Eq. ([!]), we get the following 
system of inequalities: 

e 3 -e 2 < V 12 - V 13 - V 23 , (2) 



e 2 -e 3 < V 13 - V 12 - V 23 . (3) 

We see that since V 23 > the conditions (H) and (f|) are inconsistent and, hence, the clusters 
of three LS's cannot provide double peaks. 

Let us turn to the case N = 4. Because the number of variants increases dramatically in 
this case we will restrict our search. Namely we will assume that the first LS with the lowest 
enery E\ is located in the center of an equilateral triangle while the other three LS's with 
energies e 2 < e 3 < £4 are located in the vertexes. Then the energy of repulsion takes only 
two values, V\ and V 2 (see Fig. 1). The first electron enters the system at Ep — E\ = s\. 
It is easy to see that there are only two energies competing for E 2 , which are: £\ -\- s 2 -\- V\ 
and £ 2 + £ 3 + V 2 . All the other two-electron states have higher energies. Similarly we 
conclude that there are only two candidates for E\. They are: E\ + e 2 + £3 + 2V± + V 2 and 
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£2 + ^3 + &4 + 3V 2 . Let us again assume n = 1 in the condition ([[]). Then, according to the 
general procedure, one should choose the lowest of two values for E\ and check condition 
([]]) with both candidates for E\. If we pick the first candidate for E\ and E\ + e 2 + V\ for 
El, this condition reduces to: £3 — £2 < —V 2 , which contradicts our assumption that the 
energies are ordered. Thus, the only remaining option is that from two candidates for Ef 
the second one has the lower energy. The corresponding condition for this can be written as 

e 4 -e 1 <2(V 1 -V 2 ). (4) 

Now with Ef = e 2 + £3 + £4 + 3V 2 the system of inequalities, resulting from ([[]) in a similar 
way as (H), (@), takes the form 

(e i -e 2 ) + (e 3 -e 1 )<2V 1 -3V 2 , (5) 



(e 4 - £2) - (£3 - £1) < ~V 2 . (6) 

Upon summation of these two inequalities we get 

£4 - £2 < Vi - 2V 2 . (7) 

On the other hand we have assumed that £4 > e 2 . Then the principal requirement for a 
double peak to occur reduces to: V\ > 2V 2 . If this requirement is met, the inequalities 
(H) and @ are consistent. This is illustrated in Fig. 1 where the graphical solution of 
the system (|5|), (^) is shown. In principle, one should also check that the energies of LS's, 
satisfying the system, satisfy the condition (JD as well. Note, however, that the condition 
(H) is stronger than (||), so that the latter is automatically obeyed, as can readily be seen 
from the following chain of relations: 

£ 4 - £1 < (£4 - £2) + (£3 ~ £1) < 2Vt - 3V 2 < 2(Vi - V 2 ). (8) 

Thus we have demonstrated that, within the restricted geometry considered, double 
peaks can occur provided that V\ > 2V 2 . Obviously, the relation between V\ and V 2 is 
opposite if the interaction between the localized electrons is simply the Coulomb repulsion. 
Since the distance between the first and the second LS is v^3 times smaller than the distance 
between the second and the third LS (see Fig. 1) we have V\ = V3V 2 . The situation changes 
if a metallic electrode is placed at a distance d from the plane of the localized electrons. 
Then the Coulomb interaction is modified to 

V (r) = -(-- - 1 ) . (9) 

The modified interaction falls off as 1/ r 3 and we indeed have V\ > 2V 2 as soon as the distance 
between the first and the second LS's exceeds 0.33d. 

If the system (|5]), @ is satisfied the evolution of the occupation of the cluster with 
increasing the gate voltage (Fermi level postion Ep) would be as follows. For Ep < E\ all 
four LS's are empty. At Ep = e\ the first LS in the center of the triangle gets occupied. As 
E F reaches the value E F = (e 2 + £3 + £4 + 3V 2 — £\)/2 an electron from the center moves 



4 



to one of the vertexes and two electrons arrive from the electrode and occupy two other 
vertexes. Finally at E F — e x + 3Vi the LS in the center gets occupied again. 

After realizing that double peaks are possible in principle, we turn to the question: how 
frequent are they? One could argue that double peaks are allowed only for extremely rare 
configurations that do not really occur in a finite-size sample. To answer this question we cal- 
culate the probability that in a cluster of 4 LS's the energies and distances between LS's are 
arranged in such a way that the two-electron tunneling becomes possible after the cluster is 
singly occupied. We start from the observation that the previous consideration for restricted 
geometry becomes general if, instead of (|]), we assume the model "hard core" interaction 
between the LS's: V(r) = U for r < d and V(r) = for r > d. It is important that the 
interaction takes only two values: U and zero. Then the above analysis for the equilateral 
triangle applies if the distances, r^, between the LS's satisfy the following requirements: 
r i2) r i3, r i4 < d and r 23 ,r u ,r 24: > d. If these requirements are met, the conditions for the 
double peak formation are given by Eqs.(^), (^) with V\ = U and V 2 = 0. It can be verified 
directly that for all other configurations of 4 LS's double peaks are forbidden. Then the 
calculation of the probability, V, of the occurrence of a double peak can be performed in 
a following way: we fix the position and the energy of the first LS and find the allowed 
phase volume for other three LS's. The advantage of the "hard core" interaction is that 
the intergrations over coordinates and energies are decoupled from each other. If we denote 
with g the density of LS's, then the expression for V can be presented as V = g 3 lil 2 , where 
I\ and I 2 are the phase volumes in the energy and coordinate spaces respectively: 

roo rco roc 

h= de 2 de 3 de^8(2U + e x + e 2 - e 3 - £4,)9(e 2 + £3 - £1 - £4), (10) 

J SI Js2 J £3 



h = J dr 2 J dr 3 J dr 4 6(r 23 - d)9(r 24 - d)6(r u - d)8(d - r 12 )9(d - r 13 )9(d - r 14 ). (11) 

The analytical evaluation of the first integral results in 7i = £/ 3 /3. The integral I 2 is 
obviously proportional to d 6 ; the numerical factor was found using the Monte-Carlo pro- 
cedure. Finally we obtain V = O.Gll(gUd 2 ) 3 . We see that the portion of double peaks is 
governed by the dimensionless parameter gild 2 which is the ratio of the interaction energy 
and the mean level spacing within the size of the "core." 

It is apparent that for the realistic interaction (|9|) the estimate for V emerges if one 
substitutes for U the value e 2 /nd — the Coulomb interaction at distance d. This gives 
V ~ (ge 2 d/K,) 3 . To find the numerical coefficient, the Monte-Carlo integration over the 9- 
dimensional space (6 coordinates and 3 energies) was performed using the program published 
in the book!. The program generated a random set of dimensionless (in the units of d and 
e 2 / ' nd) coordinates and energies, calculated the values E\, Ef, E\ for the interaction (§), and 
then checked condition (|I]). The numerical factor obtained is (5.1 ± 0.1) • 10~ 2 . 

The calculation of g poses a separate problem. One approach to estimate g is to assume 
that the random potential just smears the edge of the band density of states go = m/iih 2 . 
Then for energies not very deep in the tail, g is still of the order of g . The product g e 2 d/K 
can be rewritten as d/ircio, where ao = h 2 K,/me 2 is the effective Bohr radius. It may seem 
that, if d is large enough, this product could be much larger than 1 . However this is not the 
case, since with increasing d the interaction of the occupied LS's becomes important. This 
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leads to the suppression of the density of states in the vicinity of the Fermi level (Coulomb 
gap!,!). For the interaction (M) the energy dependence of g was studied both analytically^ 
and by computer simulations!^. It was shown that g(e) = 0.085(d/Ke 2 ) + ^(k 2 \e - E F \/e 4 ). 
In our problem the relevant energy scale is \e — Ep\ ~ e 2 / nd so that ge 2 d/ k is of the order 
of 1. In view of the ambiguity in g, our calculation can be considered only as an estimate 
showing that V is not small. Indeed, the only small parameter in V is the numerical factor 
0.051. This factor is, in fact, surprisingly large taking into account that it emerged as a 
result of the 9-fold integration. 

We have addressed only the energy aspects of the double peak formation. There is also 
a question about the dynamics of the process. Suppose that the occupation of a single- 
electron state requires a time r (r is inverse proportional to the probability of tunneling). 
For a two-electron state with a binding energy W it can be shown that this time increases 
dramatically and becomes of the order of r 2 W/h. The reason for the enhancement is that 
two electrons cannot tunnel sequentially because of the energy restrictions. However, if 
the temperature T is finite the sequential tunneling becomes possible due to the smearing 
of the Fermi distribution in the electrode. Assume for simplicity that both single-electron 
energies in a two-electron state are the same. Then the time required for the sequential 
occupation of the two-electron state is equal to r exp(W/2T). This time is shorter than 
r 2 W/h if T > W/[2 \n(Wr/h)]. 

Finally let us discuss the magnetic field dependence of the double peaks. Note that 
conditions ©,© (formation of a double peak) require the energies of all four LS's in the 
cluster to be rather close (roughly speaking, they should lie within the interval of the order 
of e 2 /nd). With increasing B each LS moves up in energy. Important is that the rate of 
this motion is different for different LS's. This is obvious if some of LS's originate from 
donors, located in the well, while others represent the size quantization levels in the lateral 
fluctuations of the random potential!,!. Since the rate for donors is much slower, condition 
([I]), met at B — 0, will get violated at some critical B due to the spread in the level 
positions. For higher B the double peak will split into two. Even in the case when all 
four components of the double-peak-cluster at B — are the ground states in parabolic 
confinements, their energies will depart from each other with increasing B, thus causing a 
splitting of the peak. If a confinement is characterized by the position of minimum Uq and 
the frequency of the zero-point motion ujq, the behavior of the energy level with B is given 
by s(B) = Uq + ^(^o + Co c/4) 1//2 (w c stands for the cyclotron frequency). Suppose that at 
B = two levels, £i(0) and £2(0), are anomalously close in energy (in order to participate in 
the cluster). This means that the sum Uq 1 ^ + TujOq = £i(0) is close to + TujJq = £2(0), 
while separately Uq 1 ^ and U^p can differ, say, by a factor of two. But in a strong magnetic 
field we have £\(B) —62(B) = Uq — Uq , so that departure is the typical fate of the initially 
aligned levels. 

In conclusion, we have demonstrated that double peaks in the differential capacitance 
may result from the interaction-induced correlations in the occupation numbers of LS's 
within a cluster. More conventional consequence of these correlations is that adding of one 
electron to the cluster might cause a redistribution of neighboring electrons over LS's in order 
to reduce the total energy. This process is similar to the formation of a polaron by a lattice 
surrounding an LS. Note that such a purely electronic "polaron" was studied intensively by 
Efros and Shklovskii and by Pollak and Ortuno (see e.g. the review 

sift 

in connection with 
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the density of states and inelastic transport in the Coulomb glass. Making a link to these 
works, our main result can be reformulated as follows: for interactions which fall off steeply 
enough with distance, the formation of an "electronic bipolaron" in certain compact clusters 
of LS's is energetically favorable. 

The authors are grateful to R.C Ashoori for a valuable comment. One of the authors 
(L. G.) is grateful to E.I. Rashba for illuminating discussion on the deformation fields of 
a polaron. The work at the University of Minnesota was supported by NSF Grant DMR- 
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FIGURES 

FIG. 1. Cluster of 4 LS's providing a double peak (a). Graphical solution of the system (5), 
(6). Dashed is the region within which this system is satisfied (b). 
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